Origins of Intermediate Velocity Particle Production in Heavy Ion Reactions 
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Investigation of intermediate-velocity particle production 
is performed on entrance channel mass asymmetric collisions 
of ^*Ni-|-C and ^^Ni+Au at 34.5 MeV/nucleon. Distinctions 
between prompt pre-equilibrium ejections, multiple neck rup- 
tures and an alternative phenomenon of delayed aligned asym- 
metric breakup is achieved using source reconstructed correla- 
tion observables and time-based cluster recognition in molec- 
ular dynamics simulations. 

PACS number(s): 25.70 -z, 25.70.Lm, 25.70.Mn, 
25.70.Pq 

During the last decade, reaction mechanisms studies 
of heavy ion collisions in the Fermi energy domain have 
unraveled new and interesting phenomena that can have 
major consequences on our understanding of both nuclear 
dynamics and excited nuclear matter properties. It was 
first established that binary dissipative collisions domi- 
nate the reaction cross section in this energy range for 
a large domain of entrance channel masses and asymme- 
tries Early, it has been noticed that unlike the 
lower energy case, where exclusively two primary frag- 
ments (or one fusion fragment) are formed the inter- 
mediate energy heavy ion collisions produce exotic struc- 
tures that populate the intermediate-velocity (IV) zone 
between the two main collision partners [p|-p^ . This pro- 
duction of particles and fragments at mid- velocity is usu- 
ally explained as the onset of the participant-spectator 
model of relativistic energies [T^ ]. However, in the Fermi 
energy domain, the interplay between two-body interac- 
tions and collective motions from the nuclear mean field 
modifies this simple picture and can potentially give rise 
to peculiar mechanisms leading to such an important ob- 
served IV particle production. Of those mechanisms, 
fast light particles and coalescence clusters ejected by 
nucleon-nucleon collisions and multiple neck ruptures are 
usually expected to be responsible for the main IV par- 
ticle production. 
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In this paper, we report for the first time clear distinc- 
tions between these prompt processes and an alternative 
phenomenon of delayed aligned asymmetric breakup that 
populates the intermediate-velocity zone by a deforma- 
tion rupture of mainly the heavier of the colliding part- 
ners in mass asymmetric collisions. These distinctions 
were observed experimentally with intermediate- velocity 
particle correlation analysis of Ni-f C and Ni-f Au reac- 
tions at 34.5 MeV/nucleon. 

A beam of ^^Ni accelerated at 34.5 MeV/nucleon by 
the coupled Tandem and SuperConducting Cyclotron ac- 
celerators of AECL at Chalk River bombarded alterna- 
tively a 2.4 mg/cm^ carbon target and a 2.7 mg/cm^ 
gold target. Charged particles issued from these reactions 
were detected in the CRL-Laval 47r array constituted by 
144 detectors set in ten rings concentric to the beam axis 
and covering polar angles between 3.3° and 140°. The 
first forward four rings (3.3° to 24°) are each made of 16 
plastic phoswich detectors with energy detection thresh- 
olds of 7.5 (27.5) MeV/nucleon for element identification 
of Z=l (28) particles. Between 24° and 46°, two rings 
of 16 CsI(Tl) crystals achieve isotopic resolution for Z=l 
and 2 ions and element identification for Z=3 and 4 ions 
with energy thresholds ranging from 2 to 5 MeV/nucleon. 
The Miniball forms the last four rings (46° to 87° and 93° 
to 140° ) and is constituted by PIN diode backed CsI(Tl) 
crystal detectors set in groups of 12 per ring. See |^ and 
references therein for more information on detectors and 
energy calibration. The main trigger for event recording 
was a charged particle multiplicity of at least 3 parti- 
cles. The present work is restricted to events selected 
in the off-line analysis by a total detected charge of at 
least 31 units. This ensures good reconstruction of the 
quasi-projectile and IV particles characteristics. 

In order to evaluate properties of the IV material, it 
needs to be isolated from its surrounding. This is done 
by reconstructing on an event-by-event basis the quasi- 
projectile (QP) and quasi-target (QT) emitters using a 
probabilistic reconstruction algorithm for the first nu- 
cleus and angular and velocity cuts for emissions of the 
second one. The statistical quasi-projectile reconstruc- 
tion algorithm proceeds by two steps consisting in build- 
ing probability tables for the attribution to the QP of a 
final detected particle and an application of these prob- 
abilities on an event-by-event basis. The main hypothe- 
ses on which the algorithm relies on are: (1) the event 



1 



heaviest particle is assumed to be the QP evaporation 
residue, (2) particles emitted by the QP are isotropically 
distributed around this residue, and (3) the forward emis- 
sion hemisphere of the QP is negligeably contaminated by 
emissions from other processes. Clearly, these hypothe- 
ses can only be verified in peripheral to semi-central col- 
lisions where the QP and QT are well enough separated 
in velocity space. Elements of the probability tables are 
referenced by classes of experimental impact parameters 
(bexp), residue (Zres) and emitted particle (Zp) atomic 
numbers, as well as relative velocity between them (Vrei) 
and its absolute projection along the residue velocity in 
the CM {vjli). The experimental impact parameter was 
determined using a combination of two observables that 
are both closely related to the true impact parameter 
accordin g to molecular dynamics and deep inelastic sim- 
ulations [|9[ . These observables are the total parallel mo- 
mentum of all charged particles in the forward velocity 
hemisphere of the center of mass (CM) reference frame 
(n!^™) and the anisotropy ratio of light charged particles 
(Z < 2) in that same frame (Ra) 



with a small tendancy to decrease toward the QT veloc- 
ity. 
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where Ap, At and Pp™ are respectively the projectile 
mass number, target mass number and projectile CM 
momentum, rg = 1.2/m and the constant c was experi- 
mentally adjusted to 1.2. Following hypotheses (2) and 
(3), the QP attribution probability for particles emitted 
in the forward hemisphere is fixed at unity, whereas the 
probability for the backward-emitted ones is determined 
by dividing the forward with the backward relative veloc- 
ity spectra gated by bexp, ^res, Zp and vjli- Details and 
results of this method for QP properties can be found in 
refs. |p^ , ^ , p2| . Particles detected in the Miniball were 
attributed to the QT emitter as well as particles having 
parallel velocity along the residue direction of less than 
—0.12c. All remaining final particles are assigned to the 
IV component. 

Figure ^ shows the evolution with experimental impact 
parameter of the parallel velocity distribution in the CM 
reference frame of the reaction charged products accord- 
ing to their attributed origin. In both reactions, as the 
impact parameter decreases, the mean parallel velocity 
of the QP products also decreases, reflecting the energy 
damping of the projectile. Along with this obvious ef- 
fect, particle production from the QT and IV material, 
for ^^Ni-|-Au reaction, becomes more important as the 
collision gets more violent. In the case of ^^Ni-|-C reac- 
tion, the a particle sub-structure of the ^^C target makes 
it easy to break at all detected impact parameters. Parti- 
cle production in the IV region appears also to be strong 
on the whole b^xp range. With the decrease of impact pa- 
rameter, the average velocity of the IV charged products 
in the ^*Ni-|-C system increases above uatat, whereas for 
^*Ni+Au system, it stays close to it on the whole range. 
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FIG. 1. CM Parallel velocity distributions for charged par- 
ticles of ^®Ni4-C (left) and ^**Ni-(-Au (right) reactions. The to- 
tal distributions (thick solid lines) are broken in QP (vertical 
hatches), QT (horizontal hatches) and IV (shaded histogram) 
contributions. The vertical line shows «jvjv. 

In order to get information on the time formation 
and characteristics of the IV material produced in these 
mass asymmetric heavy ion reactions, molecular dynam- 
ics simulations were performed. A semi-classical molecu- 
lar dynamics model with a spherical two-body interaction 
potential presenting realistic nucleon-nucleon cross sec- 
tions and including Coulomb interactions has been used 
to modelize semi-peripheral reactions of the ^*Ni-|-^^C 
and ^^Ni-l-^^'^Au systems at 34.5 MeV/nucleon. Focus- 
ing on the nucleon-nucleon correlations that give rise to 
the fragment formation processes, a clusterization algo- 
rithm making use of the 'minimum spanning tree in en- 
ergy space' (MSTE) criterion, has been used. Right after 
the most violent stage of the collision, an MSTE cluster 
recognition step allowed us to sort the particles as belong- 
ing to the projectile-like fragment (PLF), the target-like 
fragment (TLF), or as prompt emitted particles (PEP) 
not spatially correlated neither to the PLF nor the TLF. 
The time at which every asymptotic cluster have actually 
been emitted (ig) can be calculated within the model by 
simply tracing back the dynamical evolution of the reac- 
tion. See j23j and refs. therein for more specific details 
about the model used and the cluster time origin analysis 
of mass asymmetric reactions. 

The extracted emission times for particles originating 
from the PLF, TLF and PEP components were used to 
put into evidence the origins of different regions of the 
asymptotic velocity space. Parallel velocity spectra of 
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particles originating from these components were sepa- 
rated in figure § according to their calculated emission 
time. From panels (a)-(c), it can be seen that for the 
^^Ni-|--'^^^Au system, the IV region is populated not only 
by PEP particles but also by a supplementary contribu- 
tion coming mainly from the heavier partner of the reac- 
tion. Moreover, a different time scale can be associated 
with each kind of contribution. The PEP are mainly light 
clusters (Z < 5) that were configurationally well differen- 
tiated from the PLF and TLF very early in the evolution. 
On the other hand (see panels (e),(i),(m)), the target-like 
intermediate- velocity clusters were spatially linked to the 
TLF at PLF-TLF separation time and were emitted later 
in the evolution, mainly between 150 fm/c and 500 fm/c. 
For the ^*Ni-|--'^^C case (fourth column of fig. an im- 
portant contribution to the IV emission coming now from 
the PLF can also be recognized. Again, that contribu- 
tion is not a prompt one, but takes place in the same 
mid-range temporal interval than the one observed in the 
other reaction. 



occurring on a longer time scale. 
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FIG. 2. Histograms of asymptotic distributions. In 
the first, second and third columns, the contributions of the 
TLF, PEP, and PLF emissions have been selectively shadowed 
for the '^*Ni-|-^^''Au reaction (7fm < 6 < 8fm). The fourth 
column shows the respective calculation for the ^^Ni-I-^'^C re- 
action (3fm < b < 4.5fm) shadowing the PLF emissions. Figs. 
(a)-(d) show total distributions. In Figs, (e)-(h), (i)-(l), and 
(m)-(p) the contribution of particles with emission times of 
U < 150 fm/c, 150 fm/c< < 500 fm/c, and te > 500 fm/c 
after the collision are shown respectively. 

These findings suggest two scenarios for the origin of 
IV particles. On one hand, the prompt emission at the 
highly coUisional stage of the reaction of mainly PEP. 
On the other, a mechanism where a dynamical defor- 
mation of the heavy partner of the reaction develops, 
eventually leading to a scission followed by a Coulom- 
bian push that projects the emitted particles towards 
the intermediate-velocity range, this second alternative 




' Ni + Au at 34.5 MeV/nueleon 



6 fm< b < 9 fm 
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FIG. 3. Oaxis vs Vrei logarithmic contour plots for IV 
attributed particles in the ^**Ni-|-C (left) and ^*Ni+Au (right) 
systems for two specific ranges of bexp- 

To investigate these possible origins of IV particles in 
the data, a study has been done of the angular relations 
between their emission direction (evaluated from their ve- 
locity vector in the reconstructed center of mass of all IV 
particles in the event) and the reaction axis, defined by 
the direction of the velocity vector of the reconstructed 
QP in the system CM reference frame. Figure ^ shows 
the distribution of the angle between those two vectors 
as a function of the relative velocity between the IV par- 
ticle and the QP residue, for particles attributed to the 
IV component. Hydrogen, Helium and IMF {Z = 3 — 7) 
contributions are shown on the figure. From collisions 
of the ^'^Ni-I-Au system, it is possible to notice firstly 
the nearly forward-backward symmetry of the distribu- 
tions around 9axis = 90° for hydrogens up to IMF's, as 
well as a sizeable focusing oi Z — 1 particles at 90°. 
This can be understood as a fast emission of particles 
and excited clusters from the overlapping nuclear matter. 
Excited clusters will decay on a variable time scale and 
most of the final products will be released in an almost 
field free region, therefore with isotropic angular distri- 
bution. However, the prompt ejected particles following 
early nucleon-nucleon collisions in the overlap region will 
be subject to the Coulomb field of the two heavy nuclei 
as well as to blocking from the saturated quantum phase 
space regions. These particles will therefore be predomi- 
nantly emitted on a plane perpendicular to the reaction 
axis. In the ^^Ni-f C reaction, forward-backward sym- 
metry around daxis = 90° is destroyed, as can be seen 
from the left panels of figure ^. Along with a narrower 
distribution, an enhancement of forward directed IV par- 
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tides with small relative velocity with the QP residue is 
observed for Z = 2 particles and IMF's. This could be 
a signature for the presence of a supplementary process 
that populates the IV region from an asymmetric and 
aligned breakup of the QP. This production process of IV 
particles seems to be less dominant in the ^*Ni+Au re- 
action since no such enhancement is observed. These ex- 
perimental results go along the molecular dynamics pre- 
dictions that aligned asymmetric breakup occurs mainly 
on the heavier partner side of the IV region ||2^. The 
fact that the effect is not experimentally seen on the QT 
side of the ^^Ni-|-Au system is just a result of inadequate 
detection thresholds for particles in this region. 
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FIG. 4. a — a correlation functions broken in longitudinal 
(solid line) and transverse (dashed line) contributions for IV 
attributed particles in the ^**Ni-(-C (left) and ^**Ni+Au (right) 
systems. 

Origins of IV particles have been studied further by 
two-particle correlation functions with directional cuts. 
A similar method was originally used for proton-proton 
[||,|5| and IMF-IMF Q correlation functions in order 
to probe the space-time extent of a compound source 
in central collisions. In our case, a particles have been 
chosen for this study for their abundance in the IV region 
and for the reason that the first unstable excited state of 
the two-a particles decaying ^Be nucleus has a lifetime of 
about 130 fm/c (F « 1.5 MeV), which is about the time 
it takes to the aligned asymmetric breakup process to 
take place. The experimental a — a correlation function 
is defined as 

l + (2) 

jyback[q) 

where Ncomcil) is the coincident a pairs yield with rela- 
tive momentum q and where Nbackiq) is the background 



yield constructed by means of a modified event-mixing 
technique for IV attributed particles. Since final par- 
ticles in each event have their velocity vectors related 
to the reaction axis, a strict application of the standard 
event-mixing technique introduces unreal contributions 
in the background yield for pair of events with different 
reaction axis directions. Therefore, by taking advantage 
of the cylindrical symmetry of the detector apparatus, 
an improved event-mixing technique was used which in- 
volved a rotation of the second event of a background pair 
in the plane perpendicular to the beam, so as to make 
coincide the (f> angle of the two event reaction axis. Cor- 
relation functions constructed with this technique were 
measured for two directional cuts in the emission direc- 
tion of the two a particles relative to the momentum 
direction of the pair in the QP reference frame. Lon- 
gitudinal correlation functions stand for angles between 
the two directions ranging from 0° up to 60°, whereas 
transverse correlation functions refer to angles from 60° 
to 90°. Figure ^ shows these correlation functions for 
the two studied systems and for different impact param- 
eters. From the figure, the ground state and the first ex- 
cited state of the particle unstable *Be nucleus are clearly 
seen at relative momenta of 25 and 100 MeV/c respec- 
tively. The weaker transverse correlation function for the 
^Be ground state at low bexp was expected from simple 
Coulomb trajectory calculations including *Be decay and 
average experimental velocity space relations for QP and 
IV particles. It is due to the small relative angle between 
the emitted a particles in the laboratory reference frame 
and to the non-negligeable size of the detectors leading to 
double hits that are discriminated out of the correlation 
function analysis. However, the calculation showed that 
multiple hits contribute to a negligeable amount in the 
difference between longitudinal and transverse breakup 
for the first excited state decay in the IV region. 

The approximately equal strengths observed for the 
longitudinal and transverse correlation functions in the 
first excited state region of interest (80 < q < 120MeV/c) 
of the ^^Ni-|-Au reaction emphasize that ^Be resonances 
in the IV region bears no alignement relations with the 
QP and are therefore probably isolated from it in space- 
time at the point of their breakup. This can be un- 
derstood as the fast QP leaves the interaction region 
very early without dragging nuclear material that can 
be pushed back in the IV region later on. However, 
differences between the strengths of first excited state 
transverse and longitudinal correlation functions are ob- 
served in the ^^Ni-|-C reaction for collisions of experimen- 
tal impact parameter below 6fm. The Coulomb trajec- 
tory calculation of a system, composed initially of a ^^Ni 
QP emitting from its surface a ^Be nucleus in its first 
excited state, has shown that tidal forces were strong 
enough to flip about 50% of initially emitted longitudi- 
nal a pairs in final transverse ones, whereas only 5% of 
the reverse case occur. This effect expected for isotropic 
a particle decay from ^Be in the field of the remain- 
ing QP is however not observed experimentally. On the 



4 



contrary, longitudinal decay appears to be dominant in 
mid-peripheral reactions. Therefore, this experimentally 
observed stronger longitudinal correlation function tends 
to imply an important initial polarization in the emission 
of first excited state ^Be resonances in the longitudinal 
direction. This can possibly occur from a deformation of 
the QP in the reaction axis direction, thus constraining 
alignement conditions for the ^Be resonance formation. 
In mid-peripheral reactions, these more important defor- 
mations will lead the QP to a breakup configuration that 
will ultimately push back some of the material in the IV 
region. 

This deformation breakup effect populating the IV re- 
gion with light fragment and particles is associated to 
the delayed aligned asymmetric breakup observed in the 
molecular dynamics simulation p3| and in the experi- 
mental Oaxis distributions for the Ni-t-C system. More 
important deformations are therefore expected on the 
heavier partner side in a mass asymmetric heavy ion reac- 
tion. This can be a result of the different surface bound- 
ary conditions at the two poles of the reseparating sys- 
tem. An asymmetrically shaped neck will breakup earlier 
on the lighter partner side and will stay attached to the 
heavier partner for a longer time, potentially leading it 
to breakup in a delayed aligned asymmetric fashion. 

In conclusion, we have demonstrated in this paper the 
existence of two mainly different origins for intermediate- 
velocity particle and fragment production. The use of a 
source reconstruction technique has permitted to study 
on an event-by-event basis the characteristic angular 
emission pattern of the IV particles as well as their rela- 
tion with the remaining QP in two-particle correlations. 
Important differences between two entrance channel mass 
asymmetric reactions were observed on the QP side of 
the intermediate- velocity space. With help of time-based 
cluster recognition algorithm applied to molecular dy- 
namics simulation p3| , it has been possible to determine 
the time scales associated to two different phenomena. 
The first origin is related to prompt nucleon-nucleon col- 
lisions that occur in the overlap zone of the two colliding 
nuclei. These processes will eject light particles and ex- 
cited clusters out of the overlap on a very short time 
scale of the order of the reseparation time. Excited clus- 
ters ejected at this stage will however emit particles on a 
longer time scale. The second origin of IV particle pro- 
duction is related to the collective motion of nucleons at 
the perturbed ends of the QP and QT. Important de- 
formations will be carried by the heavier partner of the 
collision and will lead it to a mass asymmetric breakup 
aligned along the reseparation axis. This is expected to 
happen after a delay of the order of 150-500 fm/c. In en- 
trance channel mass symmetric collisions, molecular dy- 
namics calculations predict this second origin of IV parti- 
cles to happen on both QP and QT. Recent results for the 
heavier system -'^^^Sn-t-^'^Nb at 29.5 MeV/nucleon pre- 
senting the necessity of adding a surface emission compo- 
nent to the neck component at mid-velocity are compat- 
ible with this expectation M . It can also be noted that 



the phenomenon of delayed aligned asymmetric breakup 
is potentially related to previously observed processes of 
dynamical projectile splitting ]29| ] and fast asymmetric 
fission |3C| ] at lower bombarding energies. It will be in- 
teresting in future analysis and experiments to study the 
isospin dependance of the IV products according to their 
specific origin, in order to disentangle shape, thermal 
and chemical equilibration processes. In fact, a possible 
chemical motion in reaction to the slow surface/volume 
deformations involved in the delayed aligned asymmetric 
breakup mechanism could explain in part the observa- 
tion of higher N/Z ratio of mid-rapidity IMF's relative 
to total systems with no initial neutron enrichment pSl . 
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